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ABSTRACT: Due to the rapid growth of electricity demand and the 
absence of regulations, sulfur dioxide (S0 2 ) emissions from coal-fired 
power plants in India have increased notably in the past decade. In this 
study, we present the first interannual comparison of S0 2 emissions and 
the satellite S0 2 observations from the Ozone Monitoring Instrument 
(OMl) for Indian coal-fired power plants during the OMI era of 2005— 

2012. A detailed unit-based inventory is developed for the Indian coal- 
fired power sector, and results show that its S0 2 emissions increased 
dramatically by 71% during 2005—2012. Using the oversampling 
technique, yearly high-resolution OMI maps for the whole domain of 
India are created, and they reveal a continuous increase in S0 2 columns 
over India. Power plant regions with annual S0 2 emissions greater than 50 
Gg year -1 produce statistically significant OMI signals, and a high 
correlation ( R = 0.93) is found between S0 2 emissions and OMI-observed S0 2 burdens. Contrary to the decreasing trend of 
national mean S0 2 concentrations reported by the Indian Government, both the total OMI-observed S0 2 and annual average 
S0 2 concentrations in coal-fired power plant regions increased by >60% during 2005—2012, implying the air quality monitoring 
network needs to be optimized to reflect the true S0 2 situation in India. 
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■ INTRODUCTION 

Sulfur dioxide (S0 2 ) is a major air pollutant that enters the 
atmosphere through natural (e.g., volcanic eruptions, wildfires) 
and anthropogenic (e.g., combustion and release of sulfur- 
containing fuels and materials) processes. S0 2 and its 
atmospheric products (e.g., sulfate, sulfuric acid) are 
detrimental to human health, harmful to ecosystems, and 
responsible for several environmental problems such as acid 
deposition, smog formation, and visibility degradation. As fine 
particles, sulfate aerosols have major impacts on climate 
because they can reflect incoming solar radiation and act as 
cloud condensation nuclei (CCN) that modify the micro- 
physical properties of clouds. S0 2 is also a regional and global 
issue because sulfate has a long atmospheric lifetime and can be 
transported long distances. Driven by the rapid economic 
development, S0 2 emissions from India have been continu- 
ously increasing over the past two decades. 1 7 On the basis of 
the most recent emission estimates, 4 ’ 8 ' 9 India has surpassed the 
U.S. to be the world’s second largest S0 2 emitting country, 
after China, since 2010. The coal-fired power sector is the 
single largest contributor, accounting for ~50% of the national 
S0 2 emissions and ~70% of the emission increment during 
1996— 2010. 4 For power plants, an accurate emission inventory 
at the plant/unit level is critical for atmospheric chemical 


transport models, as it can improve the model simulations 
significantly. 10,11 Although there have been some published 
studies reporting S0 2 emissions from the Indian coal-fired 
power sector, few of them used actual activity rates and 
emission factors at the plant/unit level 1 -3,s- ' or presented the 
year-by-year emission trend for the most recent years. 4,12-14 
Following our previous work on NO* emissions in India, 15 we 
use a unit-based methodology to develop the S0 2 inventory for 
the Indian coal-fired power sector in this study. 

Compared to primary carbonaceous aerosols (i.e., black 
carbon and organic carbon), S0 2 problems are not given 
sufficient attention by the Indian Government. Currently, there 
is no SO, emission control regulation for coal-fired power 
plants because of the relatively low sulfur content of coals. On 
the basis of Indian government official reports, there was only 
one out of 361 monitoring stations recording an annual mean 
S0 2 concentration higher than the national air quality standard 
of 50 /tg m -3 in 2010, and there has been a decline in S0 2 levels 
since 2000. 16,17 This is inconsistent with the fast increase in 
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Indian S0 2 emissions, and an independent data source is 
needed to verify the real S0 2 situation in India. Now, satellite 
remote sensing instruments such as the Ozone Monitoring 
Instrument (OMl) make routine S0 2 observations at high 
temporal and spatial resolution, providing valuable information 
about surface S0 2 sources. 18 Compared to nitrogen dioxide 
(N0 2 ) retrievals, satellite S0 2 signals are very noisy. In the past 
few years, satellite S0 2 instruments have been only used to 
detect the global distributions, 19 ’ 20 plumes, 21 ' 22 and trends 23,24 
of very large anthropogenic S0 2 sources. Oversampling smaller 
pixels from the center of the swath to fine grids (~3 km wide), 
de Foy et al. 22 and Russell et al. 2S identified plume features that 
were not visible on coarser grids. Averaging a large number of 
OMI pixels around the sources over a period of three years, 
Fioletov et al. recently identified the S0 2 plumes 26 and 
quantified the S0 2 emissions from individual U.S. power 
plants 27 that were typically undetectable using standard analysis 
techniques. Using the same method, S0 2 emissions from 
Canadian oil sands were also studied. 28,29 In this work, we 
apply a similar oversampling technique but with improved OMI 
retrievals to India and, for the first time, study the relationship 
between OMI S0 2 observations and S0 2 emissions over Indian 
coal-fired power plants during the OMI era of 2005—2012. We 
also show that OMI has the capability to monitor the 
interannual trend of S0 2 emissions from these point sources. 

■ METHODS AND DATA SETS 

Unit-Based S0 2 Emission Inventory. In our previous 
work, we developed a bottom-up, unit-based inventory for NO* 
emissions from Indian public coal-fired power plants during 
1996— 2010. 15 In this study, we adapt the system to calculate 
annual S0 2 emissions and extend the inventory to the year 
2012. In particular, we include a large number of captive 
(privately owned) coal-fired power plants that were not taken 
into account in the previous NO* study. In total, there are 165 
coal-fired power plants containing more than 720 units in the 
updated database. On the basis of the data availability, unit-level 
or plant-level information was collected, including geographical 
location, boiler size (i.e., capacity) and type, coal type (i.e., hard 
coal or lignite) and sulfur content, electricity generation, 
specific coal consumption, S0 2 control technology, the exact 
time when the unit came into operation and/or retired, etc. 
Total S0 2 emission (E) from coal-fired power plants for year i 
is estimated by the following equation: 

E i = Z Z Z 2-G ij; . fc -SCC i ijik -S j!k il - SR 7 .)-(1 - 17 ,) 

j k i 

( 1 ) 

where j, k, and 1 represent coal type, generation unit, and 
emission control technology, respectively; G is the electricity 
generation; SCC is the specific coal consumption per unit 
electricity generation; S is the sulfur content; SR is the sulfur 
retention in ash; and t'l is the removal efficiency of S0 2 control 
technology. 

For public coal-fired power plants, most of the unit-level or 
plant-level activity rates were derived from a series of the 
Performance Review of Thermal Power Stations 30 published by 
the Central Electricity Authority (CEA), Ministry of Power of 
India. For captive plants, such information was collected 
through various annual reports of industrial companies 
operating these plants. The exact locations of power plants 
were obtained from the Global Energy Observatory (GEO, 
http://globalenergyobservatory.org/index.php) and/or verified 


directly through Google Earth. As shown in Figure SI of the 
Supporting Information (SI), the electricity generation and coal 
consumption in Indian coal-fired power plants increased 
dramatically by 61% and 64%, respectively, from 2005 to 
2012. The predominant sizes of generating units are between 
200 and 300 MW, and these units contribute about 50% of 
both the electricity generation and the coal consumption. In 
terms of the incremental growth, however, units with capacity 
>300 MW are the biggest contributors, accounting for 70% of 
the electricity generation growth and 63% of the coal 
consumption growth. The activity rates derived from this 
study match the official statistics from the Government of 
India 30 and the International Energy Agency (IEA) 31 very well 
(differences <5%), further indicating that the current database 
covers nearly all public and most captive coal-fired generating 
units in India. 

Plant-specific coal sulfur contents ( S ) were compiled from a 
large number of reports, including Mittal et al., GEO, IEA 
Coal Power Database, 32 and various power project reports. 
Where plant-specific S values were not available, the state- 
average S values reported by Reddy and Venkataraman 14 were 
used. The national average S value of coals (including lignite) 
was calculated to be 0.55% in 2005, increasing to 0.58% in 
2012, consistent with previous studies at the national level. 4 ’ 6 
Using the default values recommended by the IPCC, 33 the SR 
ratios were assumed to be 5% for hard-coal-fired power units 
and 30% for lignite-fired ones. The effects of combustion 
technology and boiler age on SR ratios were not taken into 
account since relevant information is not available. Currently, 
S0 2 emissions are not regulated in India for coal-fired power 
plants. The Indian government only stipulates that space must 
be provided for installation of a flue-gas desulfurization (FGD) 
device for new generating units with capacity 500 MW and 
above, in case stringent regulations are implemented in the 
future. 34,35 Because the installation and operation of FGD is 
costly, it is understandable that the Indian coal-fired power 
plants would only install FGD devices if and when they are 
mandatory. In this work, we assume the S0 2 emissions are 
uncontrolled for all plants except the Trombay and Dahanu 
power plants, where seawater-based FGD devices (>/ = 90%) 
were reported to be operating. 36 The capacity of these two 
plants accounts for ~2% of the total coal-based capacity in 
India during 2005—2012. 

OMI S0 2 Retrievals and Processing. The OMI is a nadir- 
viewing imaging charge-coupled device spectrometer that has 
been measuring the earthshine radiance and the solar irradiance 
in the ultraviolet— visible range and providing aerosol and trace- 
gas observations, including S0 2 , with daily global coverage 
since October 2004. 36 It is aboard the National Aeronautics and 
Space Administration (NASA)’s EOS/Aura satellite, which is in 
a sun-synchronous ascending polar orbit at 705 km altitude 
with 13:45 local equator-crossing time. In this study, we use the 
planetary boundary layer (PBL) S0 2 data in the OMS02 Level- 
2 product, developed specifically for anthropogenic S0 2 
pollution sources. 3 ,38 Daily pixel retrievals were filtered to 
remove data with large solar zenith angle (>70 degrees) or 
relatively high radiative cloud fraction (>0.2) and terrain height 
(>2 km). Pixels at swath edges (10 pixels on each side) were 
excluded to limit the across-track pixel width to ~40 km. Since 
June 2007, some cross-track positions were affected by row 
anomalies (RA) due to the partial external blockage of the 
radiance port on the OMI instrument (http://www.knmi.nl/ 
omi/research/product/rowanomaly-background.php). We dy- 
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namically removed these affected pixels based on the RA flags 
provided in the 0MS02 Level-2 data. We also rejected pixels 
with S0 2 observations higher than 5 Dobson Unit (DU, 2.69 X 
10 16 molecules cm -2 ) to exclude cases of transient volcanic 

so 2 . 27 

For valid pixels, a series of further corrections were applied to 
the originally retrieved PBL S0 2 values. First, OMI PBL S0 2 
vertical columns were converted back to slant columns by 
multiplying by the predefined air mass factor (AMF) of 0.36. 
Then, the latitude-dependent Pacific average slant columns 
were subtracted from all valid pixels on a daily basis to remove 
the negative offset and the sharp gradient at 30°N that is caused 
by the switch of the a priori ozone profile at ±30° latitude 
during the retrieval process (the so-called “Pacific sector 
correction”). 20 The corrected slant columns were then further 
divided by the GEOS-Chem model precomputed clear-sky 
monthly local AMF to obtain the S0 2 vertical columns. 20 Local 
AMF correction, Pacific sector correction, and dynamic 
correction for row anomaly are the three most critical steps 
in developing the gridded OMI S0 2 Level-3 product (i.e., 
0MS02e at a resolution of 0.25 degree), and they significantly 
improve the quality of OMI S0 2 retrievals. 20 ' 39 

The smallest size of the OMI pixel is 13 X 24 km 2 , which is 
too large to observe the S0 2 spatial distribution near a point 
source. Averaging a large number of individual pixels centered 
within a several km radius from the source, Fioletov et al. 27 
found that statistically significant signals can be obtained for 
sources with annual S0 2 emissions greater than ~70 Gg year -1 . 
Using the same approach, we oversampled the valid pixels with 
corrected vertical columns at a 2 km X 2 km grid for the whole 
domain of India to analyze the interannual variations of the S0 2 
spatial distribution over coal-fired power plants. It is reported 
that large-scale biases still remain after the above corrections 
and processing, especially in the regions with broad pollution or 
persistent volcanic degassing. 39 To remove them, local bias 
corrections were applied on a yearly basis, as recommended by 
Fioletov et al. 27 

■ RESULTS AND DISCUSSION 

S0 2 Emissions of Indian Coal-Fired Power Plants 
during 2005—2012. Figure 1 shows the annual trend of S0 2 
emissions from Indian coal-fired power plants by unit size from 
2005 to 2012. A detailed emission summary by state and by 
region is provided in Table SI of the SI. During 2005—2012, 
S0 2 emissions increased rapidly by 71% from 3354 to 5738 Gg, 
with an annual average growth rate (AAGR) of 8.0%. This 
dramatic change reflects rapid economic and social develop- 
ment and no S0 2 emission control in the Indian power 
sector. 4,13 In terms of boiler size, coal-fired units with capacity 
>200 MW are the major contributors, accounting for ~73% of 
the total emissions. In particular, S0 2 emissions from units with 
capacity larger than 300 MW increased by 1327 Gg from 2005 
to 2012, contributing 56% to the emission increment. Since the 
current OMI S0 2 products can only capture strong 
anthropogenic point or area plumes, 18 the fact that the 
majority of S0 2 emission increment in India is from newly 
installed large power units provides a favorable condition for 
observing the interannual increasing trend of Indian S0 2 
emissions with OMI retrievals, and this will be discussed in 
the next section. 

Geographically, the distribution of coal-fired power plants is 
relatively uniform across the country except for the northeast- 
ern states (Figure 2a). Uttar Pradesh, Orissa, Gujarat, 
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Figure 1 . S0 2 emission estimates of Indian coal-fired power plants by 
unit size (bars) during 2005—2012 and their comparison with previous 
studies. Emissions of INTEX-B, 7 Garg et al., 1 and EDGAR4.2 2 are 
estimates for the whole thermal power sector (i.e., including coal-, oil-, 
and gas-fired power plants). 


Chhattisgarh, Maharashtra, and Tamil Nadu are the six states 
with power-plant emissions >500 Gg S0 2 in 2012, and 
together, they accounted for 60% of total emissions in that 
year. The emission increment is highest in Chhattisgarh (330 
Gg in eight years), followed by Gujarat (318 Gg) and Orissa 
(300 Gg). Delhi is the only state where the S0 2 emissions in 
the power sector have decreased in the past few years, mainly 
due to the conversion of coal-fired plants to gas-fired ones and 
the relocation of big power plants outside of the city of New 
Delhi. 

The uncertainties of emission estimates were determined by 
a Monte Carlo approach with the Crystal Ball software, and 10 
000 simulations were performed, as in previous work. 4,15 We 
applied normal distributions with uncertainties (hereinafter, 
95% confidence intervals) of 5% and 10% to coal consumption 
data and sulfur contents, respectively. 4 ' 15 For SR ratios and the 
FGD S0 2 removal efficiency, uniform distributions were 
assumed in the range of minimum and maximum values 
reported in the literature. 4 ' 6,14 ' 40 Results show that the 
uncertainties of estimated S0 2 emissions are ~±13% around 
the mean. They are lower than in our previous work 4 ' 40 because 
of the detailed unit-based methodology and reliable unit/plant- 
level information used. Figure 1 also compares the estimated 
emissions to other inventories. 1-7,13 Clearly, all of the estimates 
show increasing trends. The AAGR of our emissions during 
2005—2010 is 6.7%, which is in line with values of 6.2%— 7.5% 
in other inventories. 2-4,13 The current estimates are in good 
agreement with our previous ones 4 calculated from the IEA 
national statistics. The discrepancies with other studies are 
caused by various factors, including the use of higher sulfur 
content 3 or lower coal consumption, 1 ' 6 the exclusion of captive 
power plants, 13 and the inclusion of oil-fired and gas-fired 
plants 2 in their calculations. 

S0 2 Emissions of Indian Coal-Fired Power Plants 
Observed by OMI. The spatial oversampling and local bias 
correction techniques described in the methodology section 
can produce high-resolution OMI S0 2 maps, and it was 
reported that comparable S0 2 sources have similar mean OMI 
S0 2 values. 27 This implies that top-down OMI S0 2 
observations can be used to estimate surface emissions and 
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Figure 2. Locations of (a) Indian coal-fired power plants, (b) their average S0 2 emissions, and (a, b) annual and (c— f) seasonal averages of OMI 
S0 2 columns over India during 2005—2012. (b) Power plants in red are selected for the 2D Gaussian fit. 
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track their changes. Figure 2 displays the spatial distribution of 
mean OMI S0 2 over India from 2005 to 2012. A number of 
satellite S0 2 hot spots are observed, and they match the 
locations (Figure 2a) and the amounts of S0 2 emissions 
(Figure 2b) of large coal-fired power plants reasonably well. 
Note that there are some data artifacts at the southern foot of 
the Himalayas that may be caused by the unreliable retrievals 
over regions of extensive snow cover and high terrain height. 
Theoretically, S0 2 in the atmosphere has a chemical lifetime of 
around one day, 41 much longer than N0 2 (several hours). As 
shown in Figure 2, however, elevated S0 2 values are found only 
within ~50 km of the point sources, a distance which is smaller 
than for OMI N0 2 observations. 15 This further demonstrates 
that the current OMI S0 2 products are very noisy and only 
large point 27 or area 28 sources with intensive average emissions 
can be detected by OMI. 

Figure 2c— f also shows the spatial distribution of average 
OMI S0 2 columns by season 42 over India during 2005—2012. 
Similar to our previous findings in OMI N0 2 seasonal maps, 15 
July to September is the worst period to observe SO, from 
OMI in India, because of the insufficient satellite observations 
caused by the frequent cloud coverage and heavy rainfall in the 
monsoon season. Different from OMI NO z , however, we do 
not observe significant seasonal variations in OMI S0 2 over 
India probably because the seasonality of the S0 2 lifetime is 
weaker. In previous studies to quantify the S0 2 emissions from 
U.S. power plants 27 and Canadian oil sands, 28 ’ 29 only May- 
August S0 2 data were used as they have higher signal-to-noise 
ratio than other months due to small solar zenith angle, low 


snow cover, strong actinic flux, etc. However, the fact that India 
is in low latitudes makes all satellite measurement conditions 
favorable all year round, meaning that the all-year satellite data 
are of high quality over India. It was found in both the U.S. 27 
and Canada that averaging over 3 years (4 summer months 
per year) of data typically produces a statistically significant 
mean OMI S0 2 map at the 95% confidence level to identify 
intensive S0 2 sources. Therefore, if all-year data can be used, 
we can generate reliable yearly OMI S0 2 maps (e.g., Figure 3) 
for India instead of requiring “3-year long-term” mean maps for 
the U.S. 27 and Canada. 28 Using such maps, we can also observe 
the interannual trend of Indian S0 2 emissions, and this will be 
discussed in detail in the next section. 

To quantify the total S0 2 amount over a source, Fioletov et 
al. 27 used a two-dimensional (2D) Gaussian function to fit OMI 
S0 2 signals within a time window and radius: 


OMI s02 = a-f(x, y) 
1 


, / 1 


“ ex F 


2(1 -P 2 ) 


(x - ft) 2 


iy - 1 *,) 2 2 p{ x - /4 )0 - n y ) 


( 2 ) 


where * and y are the coordinates of the valid OMI pixel center; 
p x and fly are the coordinates of the center of the fitted 2D 
Gaussian function; and a x , a y , and p are the parameters 
determining the shape of the fitted 2D Gaussian function. Since 
the integration of f(x, y) over the whole surface equals one, the 
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Figure 3. Spatial distribution of yearly OMI S0 2 columns over India during 2005—2012. OMI S0 2 maps for years 2009—2012 are combined 
biennially to make the satellite sample amounts similar to years 2005—2008 because, due to the row anomaly, about half of the track positions were 
operational for 2009—2012 compared to 2005—2008. 


parameter a physically means the total number of S0 2 
molecules observed (or the observed S0 2 burden) near the 
source. If OMI S q 2 is in units of molecules km -2 and rr x and 
are in km, then a is in units of molecules and can be converted 
to mass units. 

In this work, we follow the same methodology. To 
quantitatively compare the S0 2 emissions from coal-fired 
power plants with OMI S0 2 vertical columns, we first select 23 
power plant regions based on the observed satellite signals. Any 
adjacent plants sharing the same S0 2 hotspot location are 
combined into a single source with total emissions set equal to 
the sum of their emissions and the location set at their 
emission-weighted center. Shown as red circles in Figure 2b, in 
total, 65 coal-fired power plants are included in this analysis, 
and they together accounted for ~69% of the total S0 2 
emissions during 2005—2012. We also calculated the 
proportion of their emissions to total S0 2 emissions within a 
radius of 50 km from weighted centers to ensure they are the 
dominant contributor to local emissions. Here, gridded 
emissions of other sources were taken from our previous 
work 4 for years 2005—2010 and scaled to 2012 based on the 
GAINS inventory. 6 Results show that coal-fired power plants 
contribute >75% to local S0 2 emissions for all 23 regions (and 
>90% for 16 regions), clearly indicating that the enhanced OMI 
S0 2 signals over these regions are primarily caused by 
emissions from coal-fired power plants. 

For each power plant area, all the valid OMI pixels within a 
60 km radius in a year were used for the 2D Gaussian fit to 
obtain the amount of OMI-observed S0 2 (i.e., a) near the 


power plant area in that year. Figure 4 shows the scatter plot of 
annual S0 2 emissions against fitted a (i.e., the annual averaged 
OMI S0 2 burden) over all selected power plant areas for the 
period 2005—2012. We found good agreement between S0 2 
emissions and OMI-observed S0 2 burden over power plant 
areas in India ( R = 0.93). Generally, power plant regions with 
annual S0 2 emissions >50 Gg year -1 produce statistically 
significant a values. This value is somewhat lower than the 
threshold of 70 Gg year -1 derived for the eastern U.S. by 
Fioletov et al., 27 probably because the additional treatments 
(including local AMF correction, Pacific correction, and 
dynamic removal of anomalous pixels) to the OMI level-2 
product improved the quality of the S0 2 retrievals. It should be 
noted that the slope of the regression line (~2.2 h) does not 
represent the chemical lifetime (typically a few days) but 
represents the effective OMI-observed dispersion time due to 
advection, deposition, and chemical conversion of S0 2 in the 
atmosphere. 27 Figure 4 also shows the fitted results obtained by 
Fioletov et al. 27 in the U.S. Since they used OMI S0 2 retrievals 
with a constant AMF of 0.36, their results cannot be compared 
with ours directly; however, we adjusted Fioletov et al.’s fitted 
line with a local AMF of about 0.7, which is averaged from the 
GEOS-Chem precomputed local AMF over India. 20 After 
correction, the slope of the fitted line for India is still smaller 
than that for the U.S., implying that the average dispersion of 
S0 2 in India is faster than in the summertime in the U.S. 
Higher temperatures, more precipitation, a more oxidizing 
atmosphere, and different wind patterns in India may all 
contribute to the shorter dispersion time. 

dx.doi.org/l 0.1 02 1 /es4039648 1 Environ. Sci. Technoi. XXXX, XXX, XXX-XXX 




Environmental Science & Technology 


Article 


S0 2 emissions (10 29 molecules hour' 1 ) 


0.5 1.1 1.6 2.1 2.7 3.2 3.8 4.3 



14.1 


11.3 Jj 
o 
E 

S 

0 

8.5 C. 
c 

01 

3 
■Q 


5.6 


2.8 


0 
CO 

1 
o 


i 


0.0 


Figure 4. Scatter plot of annual S0 2 emissions from 23 power plant 
regions against fitted a (i.e., the OMI-observed S0 2 burden) for 
2005—2012. Each point represents a yearly fitted a for a source. Actual 
OMI S0 2 measurements are used for the fitting. Error bars express the 
95% confidence intervals. 


Interannual Trends of S0 2 over Indian Coal-Fired 
Power Plants during 2005-2012. OMI S0 2 retrievals can 
be further used to study the interannual trend of S0 2 emissions 
from Indian coal-fired power plants, because one year of data is 
sufficient to generate a high-quality oversampled OMI S0 2 map 
for India. Figure 3 shows the yearly OMI S0 2 maps over India 
during 2005—2012. A clear increase in S0 2 columns can be 
seen in nearly all hotspots, and new hotspots have come into 
existence in the past few years, reflecting the fast construction 
of large power plants/units in India. The yearly sum of OMI 
S0 2 for all power plant regions as a function of the distance 
between the sources and the pixel centers is given in Figure 5. 
Although the current OMI S0 2 products can only detect the 
elevated S0 2 values within a radius of ~60 km, the sum of OMI 
signals over the hotspot centers was continuously increasing 
during 2005—2012. The interannual variations of the OMI- 
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Figure 5. Yearly sum of OMI S0 2 for all power plant regions as a 
function of the distance between the source and the pixel center. 
Actual OMI S0 2 measurements are used for the calculation. Error bars 
are the 95% confidence intervals of the sum. 


observed S0 2 burden (i.e., the sum of fitted a) and S0 2 
emissions of all selected power plant regions are shown in 
Figure 6. The total amount of S0 2 observed by the OMI 
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Figure 6. Interannual trend of S0 2 emissions from selected Indian 
coal-fired power plant regions, the OMI-observed S0 2 burden (the 
sum of fitted a and the corresponding 95% confidence intervals), 
national mean S0 2 concentrations reported by the CPCB of 
Government of India, and annual average SO, concentrations at 
selected coal-fired power plant regions. R values shown are the 
correlation coefficients with the OMI-observed S0 2 burden. 


increased by 63% from 2005 to 2012, consistent with the 
dramatic increase of S0 2 emissions in these regions ( R = 0.93). 
The high correlation coefficient further demonstrates the close 
relationship between the OMI S0 2 observations and the 
bottom-up emissions over coal-fired power plant regions. 

Figure 6 also shows the trend of national mean S0 2 
concentrations reported by the Central Pollution Control 
Board (CPCB), Government of India. 17 These data were 
derived from the CPCB’s National Air Quality Monitoring 
Programme (NAMP), a network of 320+ operating stations in 
120+ cities/towns during 2005—2010 (http://cpcbedb.nic.in). 
The national mean of S0 2 levels reported by the Indian 
Government shows a decreasing trend, which is explained by 
recent policy measures such as a reduction of the sulfur content 
of diesel, the use of cleaner fuels in metropolitan areas, a change 
in domestic fuel from coal to liquefied petroleum gas (LPG), 
etc. 16,17 However, this declining national concentration trend 
reported officially differs markedly from the trends of both 
OMI S0 2 observations and national S0 2 emissions. 2-6 ' 13 The 
reason for the difference is that most of the CPCB-NAMP 
stations are situated in large cities having few power plants. For 
further analysis, we selected those CPCB-NAMP stations that 
are located in our coal-fired power plant regions and have 
continuous S0 2 records between 2005 and 2010. In total, there 
are 70 stations in 27 cities that meet these two conditions (see 
Table S2 of the SI for the station list), and the annual average 
S0 2 concentrations for these 70 stations are shown in Figure 6. 
Clearly, the trend of surface S0 2 levels in the power plant 
regions is in excellent agreement with the OMI observations ( R 
= 0.98) and both of them increased by ~60% during 2005— 
2010. The above analysis clearly indicates that the network of 
the CPCB-NAMP stations needs to be optimized to reflect the 
true S0 2 situation across the country. Current government 
reports and statistics may give a misleading impression that the 
S0 2 situation is improving. In general, there are too many 
monitors located in city centers where S0 2 pollution might be 
reduced due to some local control measures, whereas there are 
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too few monitors in emission-increasing areas that truly 
dominate national emissions. 

We also noticed a difference between the changes in S0 2 
emissions and the changes in OMI S0 2 observations (or 
surface S0 2 concentrations) between 2005 and 2012. As shown 
in Figure 6 , the OMI S0 2 burden and the annual average S0 2 
concentrations in coal-fired power plant regions increased by 
more than 60% during 2005—2012, while S0 2 emissions from 
these selected regions increased by only 49%. This is a similar 
situation to what we previously found for N0 2 over Indian 
power plant areas after 2005: 15 the changes in surface N0 2 are 
greater than the changes in NO* emissions. For N0 2 , we 
attributed this to possible changes in the overall NO* chemistry 
over Indian power plant areas. We hypothesized that the 
continuous increase in NO* emissions enhances OH radical 
production through chemical feedback and thus decreases the 
NO, lifetime before 2005 but consumes OH radicals and thus 
increases the NO* lifetime after 2005. For SO^ it is possible 
that the N0 2 -induced decreasing OH radical concentration 
after 2005 also reduced the oxidation of S0 2 and decreased the 
conversion efficiency of S0 2 to sulfate. As a result, the S0 2 
concentrations increased at a greater rate than the emissions. 
Large industrial emitters near the power-plant regions, which 
are not taken into account in this study, may also contribute to 
the discrepancy. Further model simulations are needed to 
explore the synergistic effects of the high growth of both NO* 
and S0 2 emissions on local atmospheric chemistry, as well as 
the consequent impacts on human health, air quality, and 
regional and global climate. 
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